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ABSTRACT 

We analyze Keck HIRES observations of a Lyman limit system a± z — 2.652 toward Q2231— 00. These 
observations afford the most comprehensive study of the physical properties of a LL system to date. By 
comparing the ionic column densities for Fe+, Fe+"'", Si+, and Si^+ against calculations derived from the 
CLOUDY software package, we have strictly constrained the ionization state of this system. This has 
enabled us to calculate accurate abundances of a Lyman limit system for the first time at z > 2, e.g., 
[Fe/H] = —0.5 ± 0.1. We also derive a total hydrogen column density, log Af(H) = 20.73 ± 0.2, which is 
comparable to values observed for the damped Lya systems. 

The system is special for exhibiting CII* 1335 absorption, allowing an estimate of the electron density, 
Tie = 6.5 ± 1.3 X 10~^ cm~^. Coupling this measurement with our knowledge of the ionization state, we 
derive the following physical properties: (1) hydrogen volume density, nn = 5.9 ± 1.2 x 10~^ cm~^, (2) 
path length, = 3 ± 1.6 kpc, and (3) ionizing intensity, log J912 — —20.22 ± 0.21. We point out that a 
number of the physical properties (e.g. [Fe/H], A(H), uh) resemble those observed for the damped Lya 
systems, which suggests this system may be the photoionized analog of a damped system. The techniques 
introduced in this Letter should be applicable to a number of Lyman limit systems and therefore enable 
a survey of their chemical abundances and other physical properties. 

Subject headings: quasars:individual (Q2231— 00) — absorption lines — galaxies: abundances 



1. INTRODUCTION 



Observations of the Quasar Absorption Line (QAL) sys- 
tems provide an efficient means for probing the physical 
conditions of the early universe. Intermediate resolution 
have enabled accurate mcaaurcmcnta on the cvo 



surve 



lutioni of the number density o f the QAL eyetoms, inclu d- 
ing the damped Lya systems ( Wolfe et al. 1986|, [Wolfe et 



al. 199f: ), the Lyma n limit (LL) systems ( Sargent, Steide. 
fc Boksenberg 19*89), and the Lya forest clouds (Sargent 



et al. 1980| ). UntlTrecently, the absence of a comprehen- 



sive physical description of the various QAL systems has 
limited the impact of these surveys. With the wealth of 
data afforded by echelle spectrographs on 10m class tele- 
scopes, one is now capable of making direct measurements 
of the physical properties of individual QAL systems and, 
in turn, the role of these systems in galaxy formation and 
the evolution of the early universe. 

Over the past few years, observations of the damped Lya 
systems - QAL systems with neutral hydrogen column 
densities A^(HI) > 2 x 10^" cm~^ - have provided detailed 
measurements on their chemical composition, ionization 
state, dust content, metallicity, and kinematic character- 
istics (iPettini et al. 1994|,|Prochaska and Wolfe 1996| , |Lu et 



al. 1996| , [Prochaska fc Wolfe 1997£j ). These observations 
present vital clues to galaxy formation theory. Thus far, 
little research has focused on the physical properties of 
the Lyman limit systems, QAL systems optically thick at 
the Lyman continuum. Owing to the challenges associated 
with accurately determining the ionization state of these 
partially ionized systems, their metallicity and chemical 
abundances have been poorly constrained. Steidel (1990) 

^Visiting Astronomer, W.M. Keck Telescope. The Keck Obser- 
vatory is a joint facility of the University of California and the Cali- 
fornia Institute of Technology. 



performed a survey of the physical properties of z « 3 LL 
systems, but the majority of his observations resulted only 
in limits to the abundances and other properties. With 
Keck, thus far researchers have foc used on the kinematic 
characteristics of the Mgll systems ( Charlton fc Churchilj 



1996), the majority of which are LL systems. Their efforts 



have provided keen insight into these systems at moder- 
ate redshift z « 1, but have been somewhat limited by 
difficulties in identifying their ionization state. 

In this Letter, we present techniques aimed at deter- 
mining the physical characteristics of a single LL system, 
including measurements of the ionization state, chemical 
abundances, electron density, and the intensity of ionizing 
radiation. While the results from this single LL system 
are not necessarily representative of all LL systems, the 
techniques outlined in this letter will be applicable to fur- 
ther LL studies, particularly when A^(HI) > lO^^cm^^. 
By applying these methods to a significant number of LL 
systems at high redshift we will investigate their chemical 
history, contrast their properties with those of the damped 
Lya systems, and thereby probe physical conditions in the 
early universe. Section 2 presents the observations and 
measurements of the ionic column densities and the ion- 
ization state. In §3, we determine elemental abundances 
by making ionization corrections to the ionic column den- 
sities. We estimate the electron density, hydrogen volume 
density, and the intensity of ionizing radiation through an 
analysis of the fine structure CII* 1335 transition. Finally, 
we provide a brief discussion and summary in §4. 

2. OBSERVATIONS AND ANALYSIS 

We observed Q2231-00 (ze™ ^ 3.02) on the night of 
November 1, 1995 with HIRES ( [Vogt 1992| ) on the 10m 
W. M. Keck I Telescope for a total integration time of 2.5 
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hours. We used the C5 decker plate with 1.1" slit, standard 
2 X 1 binning on the 2048 x 2048 Tektronix CCD, and the 
kv380 filter. For reduction and calibration of the data, we 
took a 360s exposure of the standard star BD +28 4211 
and images of quartz and Th-Ar arc lamps. The data 
was reduced and wavelength calibrated with the software 
package written by Tom Barlow and continuum fit with 
the Iraf package continuum. 



ization state, metallicity, and dust depletion. Table |l| lists 
the column densities for a number of the ions over four 
velocity regions and the total integrated column densities. 
We are particularly interested in the relative ionic column 
densities of Fe+, Fe+"'", Si"*", and Si'^+ because the ratios 
Fe+/Fe++ and Si+/Si^+ constrain the ionization state of 
the system. 




log N(HI) 




lunii illllyiillli! 



Relative Velocity (km s ^) 
Fig. 1. — The Lya profile of the Lyman limit system at 2 = 2.652 
toward Q2231-00. The fit was performed with the VPFIT software 
package and corresponds to logA'^(HI) = 19.12 ± 0.15 dex which 
includes an estimate to the continuum error. 



A crucial component in the analysis of any QAL system 
is the neutral hydrogen column density, A^(HI). Figure]^ 
shows the HI Lya profile for the LL system at z = 2.652 
toward Q2231— 00. While the Hne is heavily saturated, 
the damping wings are resolved and provide a good mea- 
surement of the N{m) value, log7V(HI) = 19.12 dex. The 
curves in Figure |^ correspond to the line-profile fit (with 
estimated errors) derived from the VPFIT software pack- 
age, kindly provided by R. Carswell and J. Webb. We in- 
clude an estimated 0.15 dex error to the Af(HI) value due 
to the difficulty in determining the unabsorbed continuum 
flux in the Lya forest. 

Figure |^ presents the velocity profiles of the metal-line 
transitions observed for this LL system, with v = corre- 
sponding to z = 2.652. The dotted lines designate regions 
of the profiles contaminated by other metal lines or Lya 
forest clouds. The dashed vertical lines overplotted on the 
Fell 1608 transition indicate the velocity regions discussed 
below (Table |^). Note that the profiles track one another 






very q loaoly in velocity apace, even across significantly dif 



feren v ionization states (e.g. Fe"*^, Fo'*'^ 



and Si^^ 



s em 
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Fig. 2. — Velocity profiles of all the observed metal-line tran- 
sitions. The dotted regions indicate portions of the profile con- 
taminated by other absorption lines. The dashed vertical lines in 
Fell 1608 define the velocity regions discussed in the text. 



To identify the ionization state of this LL system system, 
we have utiliz ed the software package CLOUDY v90.04 
( Ferland 1991 ). Figure || presents the predicted ionic col- 
umn densities versus ionization parameter U f or a model 
system assum ing a Haardt-Madau spectrum (Haardt & 
Madau 199^ ), where 



phasized below, this observation is contrary t o that typ- 
ically observed for the damped Lya syst ems ( Prochaska 
and ^olie l"996| , [Wolfe fc Prochaska 19991 ). We have mea- 
sured the ionic column densities from the unsaturated and 
mildly saturated line profiles by summing the a pparent 
colu mn density {N a{v) = mecln[/c//a(w)]/7re^/A; 3avagc 
and S |3mbach 199l| ) over several velocity regions. This ap- 
proach enables one to calculate ionic column densities as 
a function of velocity to search for variations in the ion- 



U = 



0912 
CUH 



^ (2 X 10-^) -^^^^//""^ 



(1) 



with uh the volume density of Hydrogen, and ^912 and 
J912 the flux and intensity of the incident radiation at 1 
Rydberg. Note this is a modified definition of U; it was 
chosen to facilitate future comparisons with this work. 
The relative ionic column densities were calculated with 
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the CLOUDY package assuming iV(HI) = 10i^-i2cm-2, 
intrinsic solar abundances, and [Fe/H] = —0.5 dex. The 
dashed vertical lines denote the range of U values consis- 
tent with the Fe+/Fe++ ratio (0.40 ± 0.05 dex) and the 
dotted line indicates the U value corresponding to the up- 
per limit on the Si'^/Si'^+ ratio (i.e. a lower limit to U). 
These observations imply \ogU = —2.23 ± .21, which is 
the most accurate determination of ?7 in a LL system to 
date. 
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Fig. 3. — The results from a series of CLOUDY calculations 
(Ferland 1991) assuming a Haardt-Madau ionizing spectrum and 
logA'(HI) = 19.12. Plotted are the calculated mean ionic col- 
umn densities (assuming intrinsic solar abundances and [Fe/H] = 
—0.5) versus a modified ionization parameter, V , as defined in the 
text. The vertical dashed lines indicate the observed Fe+/Fe++ 
ratio ( with eiiors) and the dotted vertical line is a lower limit to 
the ioi ization parameter determined from the upper limit to the 



Si+/Si-^+ ratio 



3. PHYSICAL PROPERTIES 

In the previous section, we presented an accurate mea- 
surement of the system's ionization state by comparing 
the observed ionic column densities of H*^, Fe"*", Fe"'""'", 
Si"*", and Si^+ with the predicted results from a calcula- 
tion performed with the CLOUDY software package. We 
find the system is highly ionized with ionization fraction, 
X = H+/H = 0.97 ± 0.02. This implies that the total 
Hydrogen column density, logA'^(H) = 20.73 ±0.2, is com- 
parable to that observed in the damped Lya systems. This 
value is similar to the estimates from Steidel (1990) for his 
sample of LL systems. By making the appropriate ion- 
ization corrections from the CLOUDY results, we derive 
the elemental abundances of this system. The logarithmic 
abundance of element X rela tive to Hydrogen and nor - 
malized to solar abundances ( Anders fc Grevesse 198^ , 
[X/H] = log[N(X)/N(H)] -log[X/Hj0, are listed in the 
last column of Table 0. We find [Fe/H] = —0.5 dex which 
is higher than the typical damped Lya system and consid- 
erably higher than any other estimate to the abundances 
of a LL system at z > 2. 



To gauge the accuracy with which we have determined 
the ionization state and thereby the chemical abundances, 
we have investigated the effects of varying 7V(HI) and the 
shape of the input spectrum. We find that the 0.15 dex un- 
certainty in the iV(HI) value lends to less than a 0.1 dex 
uncertainty in [Fe/H] and the majority of other elemen- 
tal abundances, including the iV(H) value. Furthermore, 
we considered the effects of a steeper spectrum (e.g. a 
Bregman-Harrington spectrum) and find a similar vari- 
ation in the measured abundances. Examining the rela- 
tive abundances of this system, we observe little departure 
from the Solar abundance pattern. There is, however, ev- 
idence for an overabundance of Si/Fe (the lower limit to 
[Si/H] is quite conservative given the degree of saturation 
in the Si profiles) and an underabundance of Ni/Fe as ob- 
served in the majority of the damped Lya systems. In 
the damped systems, this pattern is typically interpreted 
as the result of dust depletion and/or Typell supernovae 
enrichment. Unfortunately, the SII 1253 and SII 1259 pro- 
files in this system are blended with Lya forest clouds for 
the [S/Fe] measurement is sensitive to both of these inter- 
pretations. 

The detection of the fine structure CII* 1335 transition 
allows further insight into the physical characteristics of 
this system. The fine structure line can be excited by 
several processes, but for a highly io nized system at thi s 
redshift, electron collisions dominate ( [Morris et al. 1986| ). 
Therefore, the observed iV(CII*)/A^(CII) ratio provides a 
measure of the electron density via: iV(CII*)/iV(CII) = 
3.9 X 10-2„g [1 -t- (0.22np/ne)]. Unfortunately the iV(Cir)/iV(CI 
ratio cannot be measured directly because the CII 1334 
profile is saturated over the velocity region where the CII* 
absorption is detected^. We can accurately estimate the 
7V(CII*)/A^(CII) ratio, however, provided the assumption 
that the CII 1334 and Fell 1608 profiles track one an- 
other in velocity space. Low-ion profiles al ways track one 
another in t he damped Ly a systems (e.g. Prochaska and 



Wolfe 1996| , Lu et al. 1996) and one notes all of the transi 
tions - irrespective of ionization state - trace one another 
in this system. Our approach, then, is (1) to measure 
iV(CII) /iV(Fell) in Reg4 where the effects of saturation are 
minimal and (2) correct a TV (Fell) measurement in Reg2 
by the iV(CII)/iV(FeII) ratio to estimate iV(CII) in this 
velocity region. We find 7V(CII) = 15.41 ± 0.08 cm'^ and 
correspondingly, iV(CH*)/7V(CII) = 3.09 ± 0.75 x lO'^. 
Assuming np/rie ~ 1 for this highly ionized system this 
yields an electron density, Ue = 6.5 ± 1.3 x 10~^cm~'^. 

The Uf, value coupled with our knowledge of the ioniza- 
tion state enables an estimate of the hydrogen volume den- 
sity and the intensity of the ionizing radiation. Ignoring 
the contribution to rig from metals (a reasonable assump- 
tion when X K, I) and adopting iV(He+)/7V(He) = 0.9 
and 7V(Hc++)/7V(He) = 0.1 from the CLOUDY calcu- 
lations, we find uh = 5.9 ± 1.2 x 10~^cm~'^. With 



iV(H) 



10 



20.73 



cm , this implies a path length £ 



iV(H)/n/f w 3 ± 1.6 kpc. While this length is rather un- 
certain (> 50% statistical error) and does not account for 
clumping along the sightline, it is still one of the most 
meaningful size estimates of a QAL system to date. The 

■^We expect the absorption in Reg4 (v fa 170kms~^) for the 
CII* 1335 transition is due to an intervening Lya forest cloud 
and restrict the analysis to the absorption features in Reg2 {v 
20 kms-i). 
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Table 1 



Ionic Column Densities and Abundances 



Ion 


A 


Kegl 


KegZ 


jxego 


ixeg4 


1 U 1 


fV /TTl 
[X/HJ 


HI 


1215 








19.12 ±0.2 






CTT 




OA + (^^ 






1 4 ^9 _(_ n n"? 

±4:.4:Zi zn U.UO 




^ u .ou 


CII* 


1335 




12.90 ±0.03 










01 


1302 


14.49 ±0.02 


> 15.20 


14.22 ±0.02 


13.82 ±0.02 


> 15.33 


> -0.6 


AlII 


1670 


12.35 ±0.01 


> 13.80 


12.88 ±0.02 


12.67 ±0.02 


> 13.84 


> -0.75 


Sill 


1304 


13.58 ±0.02 


> 14.87 


14.18 ±0.01 


14.03 ±0.02 


> 15.02 


> -0.55 


SilV 


1402 




> 14.50 


> 14.21 


> 14.09 


> 14.78 




Fell 


1608 


12.96 ±0.08 


14.34 ±0.01 


13.38 ±0.03 


13.43 ±0.03 


14.44 ±0.05 


-0.5 ±0.1 


Felll 


1122 


13.71 ±0.04 


14.69 ±0.03 


13.86 ±0.03 


13.83 ±0.03 


14.84 ±0.02 




Nill 


1741 










13.36 ±0.04 


-0.78 ±0.05 



"Velocity region spanning: —100 km s^^ < v < —20 km s^^ 
''Velocity region spanning: —20 km s^^ < u < 70 km s^^ 
^Velocity region spanning: 70 km s^^ < w < 140 km 
''Velocity region spanning: 140 km s^^ < v < 210 km 

■^This value was obtained with the method described in §3 



ionizing intensity, J912, is easily determined by inverting 
Equation |l|: logJ9i2 = —20.22 ± 0.21 which is higher (a 
factor o f 5 — 10) than t ypical estimates from the Proximity 
Effect ( |Lu et al. 1991 ), but not unreasonably high. Per- 
haps there is a source of local ionizing radiation (e.g. su- 
pernovae, OB stars) in addition to the extragalactic back- 
ground radiation. Finally, consider the kinematic charac- 
teristics of this LL system. First note that the width of 
the profile spans nearly 200 km s^^ . Second, we observe 
that the strongest absorption feature lies at the left edge of 
the profile and the optical depth decreases monotonically 
towards positive velocity. This edge-leading asymmetric 
trait is ch aracteristic of that observe d in the damped Lya 
systems (Prochaska & Wolfe 19971, Prochaska & Wolfe 
19980 Contrary to the damped Lya systems, however, 
the low-ion profiles trace very closely the high-ion profiles 
(Si'^"'" and Fe++). This indicates the various ionization 
states arise from the same regions within the system. In- 
terestingly, one observes a similar correspondence between 
the line profiles of multiple ionization states in the ISM 
( Savage et al. 1997 ). A future survey will determine if a 
large sample of LL systems exhibit similar traits. 

4. SUMMARY AND CONCLUSIONS 

We have presented techniques for accurately determin- 
ing the ionization state of a LL system. The technique 
requires the measurement of multiple ionization stages of 
a single element - in this case Fe+/Fe+"'" - an estimate of 
A^(HI) and a series of calculations with the CLOUDY soft- 
ware package. It should be applicable to a number of high 
redshift LL systems and will enable a survey of their ele- 
mental abundances. For this system, we find an ionization 
fraction x — 0.97 ± 0.02 and a corresponding Iron abun- 
dance, [Fe/H] — —0.5 ± 0.1 dex. The overall abundance 
pattern resembles that observed in the damped Lya sys- 
tem possibly suggesting an underlying depletion pattern 
or Typell SN enrichment. This system is unusual for ex- 
hibiting CII* 1335 absorption which provides one with an 
estimate of the electron density (rig — 6.5 x lO^^cm"^). 



Our knowledge of the ionization state allows estimates of 
the hydrogen volume density {nn = 5.9 x 10^^ cm~'^), the 
intensity of the ionizing radiation (log J912 = —20.22), and 
the path length through the system ^ = 3 ± 1.6 kpc. 

It is instructive to compare this LL system with the 
properties of the damped Lya systems. For example, we 
note the derived -/V(H) and nn values are similar to those 
observed for the damped systems. In fact, if it were not for 
the high intensity of ionizing radiation, this system would 
be primarily neutral with A^(HI) > 2 x 10^°cm~^. The 
observed abundances and kinematic characteristics further 
argue for this interpretation. It will be exciting to inves- 
tigate the connection between other LL systems and the 
damped Lya systems through a survey of LL systems at 
z > 2. This research will help develop a better under- 
standing of the connection between the QAL systems and 
galaxy formation. In addition, a direct comparison with 
numerical simulations may provide insight into the nature 
of protogalaxies in the early universe. 

We wish to acknowledge A. Wolfe, S. Buries and F. 
Hamann for helpful comments. We would also like to 
thank Bob Carswell and John Webb for kindly provid- 
ing the VPFIT package, Piero Madau for providing an 
electronic version of their UVB spectrum and G. Ferland 
for developing and distributing CLOUDY. We thank T. 
Barlow for developing and providing the HIRES data re- 
duction package. Finally, this work could not have been 
carried out without the helpful support of the Keck ob- 
serving staff. JXP was partially supported by NASA grant 
NAGW-2119 and NSF grant AST 86-9420443. 



REFERENCES 

Anders, M.V., & Grevesse, N. 1989, Geochim. Cosmo. Acta, 53, 197 

Bregman, J.N. & Harrington, J. P. 1986, ApJ, 309, 833 

Charlton, J.C. & Churchill, C.W. 1996, ApJ, 465, 631 

Ferland, G. J. 1991, Ohio State Internal Report 91-01 

Haardt, F. & Madau, P. 1996, ApJ, 461, 20 

Lu, L., Wolfe, A.M., & Turnshek, D.A. 1991, ApJ, 367, 19 



PROCHASKA 



5 



Lu, L., Sargent, W.L.W., Barlow, T.A., Churchill, C.W., & Vogt, S. 

1996b, ApJS, 107, 475 
Morris, S.L, Weymann, R.J., Foltz, C.B., Turnshek, D.A., Shectman, 

S., Price, C, & Boronson, T.A. 1986, ApJ, 310, 40 
Pcttini, M., Smith, L. J., Hunstead, R. W., and King, D. L. 1994, 

ApJ, 426, 79 

Prochaska, J. X. and Wolfe, A. M. 1996, ApJ, 470, 403 
Prochaska, J. X. & Wolfe, A. M. 1997, ApJ, 474, 140 
Prochaska, J. X. & Wolfe, A. M. 1997, ApJ, 486, 73 
Prochaska, J.X. & Wolfe, A.M. 1998, ApJ, in press 
Sargent, W.L.W., Young, P.J., Boksenberg, A., & Tytler, D. 1980, 
ApJS, 42, 41 

Sargent, W.L.W., Steidel, CO., & Boksenberg, A. 1989, ApJS, 69, 
703 

Savage, B. D. and Sembach, K. R. 1991, ApJ, 379, 245 
Savage, B.D., Sembach, K.R., & Lu, L. 1997, AJ, 113, 2158 
Steidel, C.C. 1990, ApJS, 74, 37 

Vogt, S. S. 1992, in ESQ Conf. aiid Workshop Proc. 40, 
High Resolution Spectroscopy with the VLT, ed. M.-H. Ulrich 
(Garching: ESO), p. 223 

Wolfe, A. M., Turnshek, D. A., Smith, H. E. and Cohen, R. D. 1986, 
ApJ, 61, 249 

Wolfe, A. M., Lanzetta, K. M., Foltz, C. B., and Chaffee, F. H. 1995, 

ApJ, 454, 698 
Wolfe, A.M. & Prochaska, J.X. 1999, ApJ, submitted 



